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Abstract 

I discuss the calculation of QCD jet rates in e + e~ annihilation as 
Oh | a testing ground for parton shower simulations and jet finding algo- 

jy^ ■ rithms. 

1 Introduction 

43 ' 

Ph| The production of jets of hadrons in all kinds of high-energy collisions is 

dramatic evidence of the pointlike substructure of matter. QCD predictions 
of the rates of production of different numbers of jets are well confirmed 
and provide good measurements of the fundamental coupling as- The latest 
triumph in this respect is the calculation of the 5-jet rate in e + e _ annihilation 
to next-to-leading order, i.e. 0(ceg) pQ. Figure Q] shows that calculation 
compared to data from the ALEPH experiment at LEP [2j. The observable 
shown is L 45 = — In (2/45), where y 45 is the value of the jet resolution parameter 
at which five jets are just resolved using the k t -jet algorithm[3]. There is 
good agreement over the range shown, and the uncertainty in the prediction 
is remarkably small considering this quantity is 0(a^) at leading order. The 
value of the strong coupling obtained from the NLO fit to the region L 45 < 6 
is 

X'- a s (M z )= 0.1156^, (1) 

which is in good agreement with the world average value obtained from other 
observables. 

However, looking at a wider range of 1/45 values, fig. El we see that the 
region used in the NLO fit represents only a small part of the full distribution. 
Most events have L 45 > 6, with a distribution that turns over at L 45 ~ 8, 
whereas the fixed-order prediction continues to rise more and more rapidly 
with increasing L45 (note the logarithmic vertical scale in fig. [1]). 
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Figure 1: ALEPH data |2J on the differential 5-jet rate, with the NLO pre- 
diction from ref. [I]. 
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Figure 2: ALEPH data [2] on the differential 5-jet rate, with event generator 
predictions. 
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What this means physically is that most events have a two-jet structure 
that can only be resolved into five jets by using a high-resolution jet algo- 
rithm. However, it is important to understand this internal structure of the 
jets as well as possible, for example to search for highly-boosted new particles 
whose decays might look like or be hiding inside QCD jets. To achieve better 
understanding we need progress on two fronts: 

1. Calculations of jet substructure in the region beyond the reach of fixed- 
order perturbation theory; 

2. Jet algorithms that probe jets in a way that reveals their substructure 
in informative ways. 

Although the era of LEP physics is past, e + e~ annihilation can still serve 
as a good testing ground for ideas on both these topics, as I hope to illustrate 
in the following sections. 



2 Parton showers 

The reason for the breakdown of fixed-order predictions at high L45, where 
most of the data lie, is that QCD matrix elements have soft and collinear 
singularities that give rise to logarithmic enhancement of higher-order con- 
tributions. In fact there are up to two factors of L 45 for every extra power 
of as, so if the coefficient were unity we would expect a breakdown at 
£45 ~ l/y/a~s ~ 3. As we shall see, in fact the coefficient is more like 
2/37T, which does indeed imply a breakdown at L 45 ~ 6. Ideally we would 
like to be able to sum these enhanced terms to all orders in a closed form 
that would exhibit the turnover in the distribution, as is the case for several 
other e + e~ observables. 

In ref. [3] we wrote down integral equations for generating functions that 
can be used to compute the leading and next-to-leading logarithms (NLL) in 
jet cross sections to any order. Table [1] shows the results up to These 
equations are for the jet fraction R n (y C ut), which is the fraction of events that 
have precisely n jets at resolution y cut . The differential jet rates, like the one 
in figs. [T] and [21 are obtained from them by differentiating: 



da g dR n 



ctot dy k - ltk n=k dy cut 

Thus to NLL accuracy, in the notation of table [Q 
1 da a 3 



(2) 



ycut— yk-i,k 



o~tot dy 45 2/45 



(6R 56 Ll 5 + 5R 55 L%) + 0(a 4 s ) . (3) 
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Table 1: Jet fractions in e + e — > hadrons to NLL order in L = ln(l/y cut ), 
expanded to third order in a = a^/n. 

R 2 = l + a(R 21 L + R 22 L 2 ) + a 2 (R 23 L 3 + R M L i ) + a 3 (R 25 L 5 + R 26 L 6 ) + ... 

R 2 i = 3C F /2 

R 22 = —Cp/2 

R 23 = -3C F /A-UC F C A /36 + C F N f /18 

R 2 4 = Cp/8 

^25 = 3C 3 F /lQ + llC 2 F C A /72-C 2 F N f /3Q 

R26 = -CJ./48 

R 3 = a(R 31 L + R 32 L 2 ) + a 2 (R 33 L 3 + R u L i ) + a 3 (R 35 L 5 + R 36 L 6 ) + ... 

R31 = —3C F /2 

R32 = Cf/2 

R 33 = 3C 2 F /2 + 7C F C ' A /12 - C F N f /\2 

Ru = -C%/A-C F C A /A8 

R 35 = -9CJ./16 - 137C 2 F C A /288 - 7C 2 A C F /160 + 5C F N f /72 + C F C A N f /160 

R 36 = C 3 F /16 + C 2 C A /96 + C F C 2 A /960 

R4 = o^^R^L 3 + R^L^ + a, 3 (^R^L^ + R^gL®) + . . . 

i? 43 = -3C£/4 - 5C F C A / 18 + C F N f /36 

R M = C 2 /8 + C F C A /A8 

R 45 = 9C 3 F /16 + 71C 2 F C A /lU + 217C F C 2 A /2880 -41C 2 F N f /720 -C F C A N f /120 

R 46 = -C 3 F /16 - C 2 F C A /A8 - 7CVC1/2880 

-R5 — a 3 (-R55-^ 5 + i? 56 L 6 ) + . . . 

R 55 = -3CJ./16 - 49C|C A /288 - 91(7^/2880 + llCj.N f /720 + C F C A N } /480 

i? 56 = C3 /48 + C A /96 + C F Ci/720 
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However, such fixed-order NLL predictions are not much use as they are 
invalid when L 45 is not large and need to be resummed when it is large. 
Indeed, since 6R^q = 197/270 = 0.73 while (for n/ = 5 flavours) 5i?55 = 
—7.77, the prediction dSJ) is actually negative for L 45 < 10. 

The leading double-logarithmic 'abelian' terms, i.e. those proportional to 
(aCpL 2 ) n ~ 2 , resum to an exponential form: 

RSI ~ ^ QaCVL 2 )" exp (—oCfL 2 ) (4) 

This gives the correct qualitative features of the differential distribution (j2J) 
at large L, but the numerical values are wrong, e.g. the turn-over occurs at 
£45 ~ 10. This is not surprising in view of the comparable non- abelian terms 
and large NLL corrections. 

The easiest way to resum the enhanced terms more completely is to en- 
code them in a parton shower simulation. By this I mean a sequential 1 — > 2 
parton branching process with branching probabilities of the form 

dP(a ->• be) = ^^-—P ba (z)dz (5) 
7r q 

where q is an ordered evolution variable, z measures the energy fraction in 
the branching, P ba is the corresponding DGLAP splitting function and the 
argument q' of a>s is a function of q and z in general. The integral equations of 
ref. j3] are equivalent to such a process with the following simple properties: 
the evolution variable is the angle of branching and q' is the relative transverse 
momentum. 

The HERWIG[4] event generator results shown in fig. [2] are based on a 
parton shower with precisely these properties. PYTHIA[S] also has a parton 
shower which, although organized in a different way, ought to be equivalent. 
ARIADNE [6] is based on a different approach involving colour dipoles rather 
than partons|§ All the generators correctly reproduce the main features of 
the distribution, in particular the turn-over at L 45 ~ 8. 

It should be said that the event generators include a lot of additional 
refinements, such as matching to fixed-order matrix elements at low L45 and 
modelling of hadronization. In particular the latter has quite a strong effect 
at LEP energies and introduces free parameters which can be tuned to the 
data. Nevertheless a parton shower, or equivalent, with the correct features 
is an essential component for reliable extrapolation to the higher energies 
and different processes encountered at the LHC. 

2 I should emphasise that the discussion in this paper concerning alternative evolution 
variables and the colour structure of the shower refer only to parton showers as defined 
by eq. §5§ and not to dipole showers. 
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Angular ordering is not the most convenient organization of the par- 
ton shower: physical quantities such as transverse momenta and jet masses 
have to be reconstructed from the shower variables. It would also be prefer- 
able to generate the hardest (highest transverse momentum) branchings first, 
which would make matching to fixed-order matrix elements [7] and NLO 
improvements [SI E] simpler. These considerations lead us to look at what 
happens if we order the shower in relative transverse momentum (p t ) rather 
than angle. 

Unfortunately with simple p t -ordering things start to go wrong even at 
the leading-log level as soon as gluon branching is involved. Instead of the 
results in tabled] for the LL coefficients in the 4-jet and 5-jet fractions, 

-R44 = Cp/8 + CVCa/48 , 

#56 = C 3 F /48 + C 2 C A /96 + C F C 2 A /720, (6) 

we getHI 

R& t] = C 2 /8 + C F C A /2A, 

= C% /48 + C 2 F C A /A8 + 13(7^/2880 . (7) 

We could try to fix things up by ordering in p t and rejecting branchings 
that are disordered in angle. For the 4-jet rate this cures the problem with 
gluon branching, fig.[3](c), but spoils the result for sequential quark branching, 
fig- E^b), while for the 5-jet fraction everything is wrong: 

I&' e) = 5C 2 F /A8 + C F C A /A8 , 

R^ = 7Cf/576 + 13C|Ca/1440 + C F C 2 A /960 . (8) 

To see what is going wrong, consider the (^2,^2) integration regions for 
diagrams [3](b) and (c), depicted in fig. HI Here zi,z 2 and 61,62 are the 
(smaller) gluon energy fractions and opening angles in successive branchings, 
and e = r cut . Thus in diagram E^b), p r ordering corresponds to e < 22^2 < 
Z\6\, giving the integration region A+B. However, the correct region is the 
angular-ordered one A+C. If we impose angular ordering after p r ordering, 
we get only A, i.e. a deficit in the coefficient of Cp. 

Now it happens that for this diagram the inclusion of region B compen- 
sates for the loss of C as far as the logarithms are concerned, so in this case 
Pt-ordering alone gives the same result as angular ordering. I will come back 
to this point later. 

3 Thanks to Mike Seymour for pointing out an error in my original calculation of the 
coefficient of C F C\. 
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(a) (b) (c) (d) (e) (f) 

Figure 5: Leading order diagrams for e + e~ — > 5 jets. 

In diagram Etc), p t -ordering corresponds to e < Z1Z2O2 < Z\Q\, i.e. ej z% < 
Z262 < d\i as shown on the right in fig. HI The region C has disappeared and 
the p r ordered region A+B is just too large, giving an enhanced coefficient 
of CfCa- However, because region C is not there, imposing angular ordering 
after p t -ordering is equivalent to simply angular ordering, giving the correct 
region A and hence the correct coefficient of CfCa- 

So perhaps the correct prescription for a ^-ordered shower is to angular- 
order only the g — >■ gg vertices? This corrects the 4-jet rate but in the 5-jet 
rate the coefficient of CpC\ is too small: 

j^ivugg) _ (7^/3 _|_ CVCa/48 , 

R^g) = C 3 p /4g + c 2 F C A /96 + C F C 2 A /9Q0 . (9) 

However, the reason for this is the same as before: if the gluon that branches 
a second time in fig. Mf) is the harder one coming from the first gluon 
branching, the situation is as on the left in fig. HJ and we should not angular- 
order the second gluon branching. 

In summary, the way to get the correct LL (and NLL) jet fractions, to 
all orders, from a pj-ordered parton shower is to enforce angular ordering 
with respect to the branching at which each parton was "created", where 
this means the branching at which it was the softer of the two produced [7]. 
More precisely, one should veto branchings that are disordered in angle with 
respect to their "creation". Technically, a veto means not branching but 
resetting the pt scale as if the branching had occurred. This is a common 
kind of procedure in parton shower generators anyway, for example to correct 
for flavour thresholds or higher orders in the running coupling. 

This looks like a better way to do parton shower event generation. With 
Pt-ordering one can more easily correct the prediction to NLO, or indeed to 
any fixed order in as in principle. One only has to correct the first few steps 
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(a) (b) 

Figure 6: Large- colour structure of wide-angle gluon emission associated 
with the parton branching a — > be: (a) angular-ordered shower; (b) p t -ordered 
shower. 



in the shower. Unfortunately there is a catch. Everything works fine at the 
parton level as far as the distribution in phase space is concerned, but the 
colour structure of the partonic final state is not correct. 

Coming back to fig. |U (left), we see that, compared to angular ordering, 
Prordering includes a region of softer, wide-angle gluon emission, B, in place 
of a region of harder, more collinear emission, C. What this means is that 
gluon radiation is moved around within the shower, the amount and distri- 
bution remaining the same. This is depicted schematically in fig. El where 
for simplicity we show the large- N c approximation, as used for hadronization 
in event generators. In fig. Ufa), angular ordering assigns a soft, wide-angle 
gluon, actually emitted coherently by partons b and c, to the parent parton 
a, which is reasonable because a does have the coherent sum of the colour 
charges of b and c. In contrast, p t -ordering assigns this gluon to the harder 
of b and c, in this case c, as in fig. E^b). That is reasonable as far as the 
momenta are concerned, but it spoils the colour structure by treating c as 
the colour source and neglecting the coherent contribution of b. 

The colour structure matters when one wants to interface the parton 
shower to a non-perturbative hadronization model. In the cluster model 
used by HERWIG, colour-singlet clusters are formed by splitting gluons at 
the end of the shower into qq pairs. Thus in the angular-ordered fig. |6]^a) the 
clusters connect (gb) and (be), while in p r ordered fig. Mjo) they connect (bg) 
and (gc). Similarly in the PYTHIA string hadronization model, the string 
connects a — g — b — c in fig. [H^a) but a — b — g — c in fig. Mb) . 

In conclusion, an angular-ordered parton shower sums the LL and NLL 
enhanced terms and provides partonic final states with colour structure con- 
sistent with QCD coherence. This is good for hadronization models but 
not so convenient for reconstruction of kinematics or for systematic improve- 
ment away from the soft and collinear regions. A p r ordered shower is better 
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in those respects and, with the right angular veto procedure, can give the 
correct NLL jet fractions. However the colour structure then needs to be re- 
configured according to angular ordering before the partonic final state can 
be hadronized. 

3 Jet algorithms 

Recall that the ^-algorithm for e + e~ annihilation j3] is defined in terms of 
the resolution variable 

Vij = 2mm{E?,E]}(l -cos^)/Q 2 , (10) 

where Eij are the energies of final-state objects i and j, 9ij is the angle 
between their momenta and Q is the centre-of-mass energy. The two objects 
with the smallest value of t/ij are combined into one, this is repeated until 
all Dij > y cnt , and the remaining objects are called jets. For the purpose 
of counting large logarithms of y cut , we can write this in the small- angle 
approximation 

€ij = min{Ei, E^Q^jQ > e , (11) 

where as before e = y/y cn %- 

As pointed out in ref. [ID], this is just one of a continuum of possible jet 
algorithms with resolution variable 

e ij = mm{Ef,E^}9 lJ /Q p , (12) 

where p can be any positive or negative number. In particular p — — 1 defines 
the resolution for the e + e~ analogue of the anti-fc 4 algorithm [10J, which has 
the advantage that objects are combined starting with those that have the 
highest energy rather than the lowest. 

When p < a supplementary condition is needed, otherwise infinitely 
soft emissions would be resolved. For anti-A^ we define 

eij = mm{Q/Ei ,Q/Ej} % , 

€i = eQ/E t . (13) 

Then if the smallest of the set of {e^, e^} is an we remove i from the list 
of objects to be recombined, and if < 1 we call it a jet. Otherwise we just 
throw it away. Thus every jet has an energy greater than eQ and is separated 
from other jets by an angle greater than e. The resulting LL coefficients in 
the 4- and 5- jet fractions are 

panti r~i2 In , r~i /~t la 

-K44 — ^Fl^ + ^F^A/o , 

R*f = C 3 F /6 + C 2 F C A /8 + C F C 2 A /48, (14) 
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where as before the large logarithm is defined as L = — 2 hie. We could 
introduce an angular resolution 5 different from the energy resolution e by 
multiplying ey by e/8. This would just replace In 2 e by In e In 5. 

It is easy to see that leading double-logarithmic abelian terms in the anti- 
kt jet rates resum to an exponential form with twice the exponent of the kt 
rates (@J: 

^i a + t ab) ~ i (aC F L 2 )\ W (-aC F L 2 ) (15) 

It should also be possible to resum the non-abelian and NLL terms using 
techniques like those of ref . [3] . 

4 Conclusions 

Although the era of high-energy e + e~ collider experiments is past, at least 
for a while, it is helpful to study how our tools for analysing hadronic final 
states perform in the cleaner environment of the annihilation process. 

The fct-jet algorithm has proven useful in all kinds of processes and the 
e + e~ jet rates defined in this way are a good place to test alternative re- 
summation methods, particular those involving parton showers ordered in 
different ways. We have seen that angular-ordered and pt-ordered showers 
can both be arranged to resum the leading and next-to-leading logarithms 
of the fcf-jet resolution y cut . The p r ordering option is good for matching to 
fixed-order calculations but causes some difficulties in matching to hadroniza- 
tion models at low scales, owing to its disordered colour structure. 

The rather different anti-fcj algorithm has been adopted as the preferred 
tool for jet finding at the LHC. An analogous e + e~ algorithm can be defined 
and we saw that it has a simple pattern of leading logarithms, which should 
be amenable to resummation using techniques similar to those applied to the 
k t algorithm. 
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